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Abstract 

The gauge fixed polygon model of 2+1 gravity with zero cosmological constant 
and arbitrary number of spinless point particles is reconstructed from the first order 
formalism of the theory in terms of the triad and the spin connection. The induced 
symplectic structure is calculated and shown to agree with the canonical one in 
terms of the variables. 

1 Introduction 

In the beginning of the 90's the gauge-fixed polygon model of 2 + 1 dimensional 
gravity with zero cosmological constant constructed by 't Hooft appeared in the 
literature pQ. The motivation was to show that closed timelike curves cannot ap- 
pear in a universe with the topology M 3 = M 2 x R, started off with good initial 
conditions and containing only spinless point particles. Soon it was discovered that 
the model gives rise to an interesting dynamical system with surprising complex- 
ity near the Big Crunch and that the direct geometric meaning of the phase space 
variables may enable one to make predictions about the possible appearance of the 
initial and/or final singularity 0IS]- There, the studies were based on numerous 
computer simulations, where the universes were constructed by intelligent guesses 
of the not independent phase space variables such that all constraints among them 
admit solutions. Later, it has been noticed [I], that if the Hamiltonian is taken to be 
the total curvature of the spacelike slice expressed in terms of the phase space vari- 
ables, then it generates the correct time evolution if the symplectic structure is the 
canonical one (pidqi). This was a promising discovery with respect to quantization, 
and it lead to a qualitative analysis of the possible spectrum of the variables in the 
quantum theory. However, due to the nonpolynomial form of the Hamiltonian and 
appearance of the discrete gauge artifacts of the model called "transitions" during 
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the otherwise trivial time evolution of the phase space variables, the quantization 
was not completed. 

There was more progress in quantization of pure gravity, because in jSJ E] it has 
been discovered that 2 + 1 gravity - apart from important subtleties [Zj - is a gauge 
theory of a connection, which is composed of the triad and the spin connection. The 
extension to the inclusion of particles was readily made |Hj, but the research mainly 
concentrated on the case of pure gravity (see [U] for an overview), apart from a few 
papers where path integral ^U], combinatorial quantization JTj or lattice methods 
[12] have been applied for treating particle degrees of freedom coupled to gravity. 
The smooth first order formalism was studied extensively by Unruh and Newbury 1 
[T3] (and quantized if one particle is present) by Matschull and Welling [HI HE], see 
also [16J. The idea in these approaches is a feature of three dimensions that one can 
reduce the infinite dimensional space of the physical fields to the space of a finite 
number of variables. Note also, that with a few exceptions ^TJEIj explicit results 
mainly correspond to topologies where M 2 has genus less than two or noncompact, 
even though topology change is not a priori forbidden in a quantum theory ^B] • The 
different approaches are usually difficult to compare and often they provide different 
quantizations. In particular, the relation of the 't Hooft polygon model to other 
approaches is poorly understood. It is a representation of Poincare group (ISO(2, 1)) 
structures on Minkowski space so it is naturally a second order approach to 2 + 1 
dimensional General Relativity. As opposed to the mentioned first order approaches 
[TH|llf)j. it contains an explicit time slicing in terms of piecewise flat Cauchy surfaces 
with some conical singularities which do not necessarily correspond to point particles 
(which are 3d curvature singularities). This is different both from the York gauge, 
where the slices are characterized by the value of the trace of the extrinsic curvature 
and the conformal gauge where the extrinsic curvature of the slices vanishes. The 
arising dynamical system has an explicitly written Hamiltonian for arbitrary genus 
of M 2 and arbitrary number of non static point particles. Recently, there has been 
progress in the spin foam quantization and loop quantum gravity, where spinning 
particles were also included [T[JJ I2D] • The results there are related to the polygon 
model; it can also be derived from the first order variables, and the the fundamental 
object describing space time is a decorated graph similarly to the spin networks of 
the mentioned loop and foam approaches. 

The comparison with a triangular lattice version of the first order formalism has 
been done by Waelbroeck and Zapata [21]. It is worked out for the smooth case 
in this article. The first step is a choice of graph T in space M 2 such that the the 
complement is contractible. Inside, the flat connection is given by a pure gauge, 
and the system reduces to the dynamics of the "edge" vectors and the Lorentz 
holonomies, specifying the contractible region in M 2 and the identifications at the 
boundaries. The second step is a gauge transformation to make all edge vectors lie 
in a spacelike surface. Then one has to show that only the lengths of the vectors and 

1 Thanks to P. Mcnotti for this reference. 



2 



the so called boost parameters of the holonomies are independent, the symplectic 
structure reduces to the canonical one and the remaining constraint algebra closes, 
no second class constraints appear. 

In the next section an introduction is given to the model with special emphasis 
on the particle content and the phase space structure. Section 3. is a brief summary 
of the first order formalism, Section 4. contains the reduction of fundamental fields 
to a finite number of covariant variables, which were essentially done in ^B] for 
trivial tangent bundle of the spacelike slice. Finally the last section contains the 
main result of the paper: the reduction of the covariant variables to those of the 't 
Hooft polygon model, the induced symplectic structure is recovered both directly 
from the fields and from the covariant variables and the dynamics is described. 

2 Polygon model 

In this section an introduction is given to the model based on the so called one- 
polygon tessellation. It applies to space times which admit a slicing in terms of 
one spacelike Euclidean polygon. In other words the space M 2 has a piecewise 
flat 2-metric with some conical singularities. There exists no rigorous proof that 
all solutions to the Einstein equations with pointlike spinless sources admit such 
slicing. However, there is ongoing work establishing this result in the mathematics 
literature [21] for the case of matter free universes, see also [22] • In the general case 
one needs more polygons with appropriate pasting conditions. However, the phase 
space structure of the model can be described in a compact way with one polygon, 
that is with a graph embedded in M 2 such that its complement has one component. 
The discussion of the general case is technically more difficult, but it contains no 
new physics. The phase space variables are associated to the edges of this graph 
and the edges of the dual graph, which are closed curves in this case. 

One may need to consult the original papers |2| for throrough physical inter- 
pretation and intuition behind the construction. Here the introduction is similar to 
reference |2Sj, where only the matter free case was discussed, although the emphasis 
here is on the the gauge fixing: the fact that the spacelike slice is a planar poly- 
gon makes all but one scalar parameter redundant of both the edge vector and the 
Lorentz holonomy. 

Suppose that the M 3 = S 9 at x / is a globally hyperbolic flat Lorentzian manifold 
with a Cauchy surface of the topology £ 9 ,at: a Riemann surface of genus g with N 
number of punctures corresponding to (spinless) particles and / is an interval. For 
g = suppose that N > 2 and for g — 1 suppose that N > 0. In the following, we 
describe the initial value Cauchy surface and its parametrization, then explain how 
it arises. 
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2.1 Phase space 



The phase space variables of the model are encoded in a decorated fat graph T on 
the surface £ 9 ,iv (to be called E from now on). The graph satisfies two properties. 

• The complement P = T\T is simply connected and has only one component. 

• It is a trivalent graph except for N edges ending at the punctures. 

In other words V has one face F — 1, then consequently the number of its edges is 
E = 6g — 3 + 2N and the number of its vertices is V = 4g — 2 + 2N. We denote 
the vertices to which three edges are incident by 3-vertices and those ending at 
punctures by 1-vertices. Then for N = one has 2E = 3V and F + V — E = 2 — 2g. 
The induced metric on P is Euclidean, it is a geodesic polygon. The initial value 
surface is P modulo gluing prescribed by the (cut) fat graph which is the boundary: 
dP = T. Examples are given in Fig^ In other words the surface S is endowed with 
a locally flat metric, with conical singularities at the vertices of T. The decoration of 
r is an assignment of one real and one positive number 77 and L, respectively, to each 
edge which characterizes the Euclidean structure completely in the following way: 
The numbers L denote the lengths of the edges. The numbers rj encode the angles 
of the polygon in such a way that any three angles corresponding to a 3-vertex of 
T are (essentially) the angles of a hyperbolic triangle with edges 2rji, 2r)j , 2rj k , where 
i, j, k denote the edges incident to the vertex. The angles at 1-vertices of T are 
determined by the mass of a corresponding particle and the parameter 77 associated 
to the corresponding edge. The precise relations will be specified below. We will 
need the dual graph 7 of T. It has one vertex p inside P and each edge d is the 
simple closed curve on £ from p to p such that it intersects edge i of T and it 
intersects only that edge. We choose an orientation for all dual edges of 7 and 
extend the decoration of T by denoting the side of edge i where the dual edge is 
oriented toward the other side by i + and the other side by That is, all edges 
of the 12g — 6 + 4N sided polygon P has a label and those with i + and i- are 
identified with opposite orientation. The quotient is the manifold E. Remember, 
that a 1-vertex is a puncture, that is the dual edges corresponding to the edges of 
T ending at one vertices are also noncontractible. Then the homotopy class (with 
endpoints fixed) of a dual edge is unique, no two distinct dual edges are in the same 
homotopy class and no dual edge is in the trivial class. E\7 consists of N discs and 
4g — 2 + iV triangles. 

Where does the above structure come from? A solution to the Einstein equations 
is a locally flat manifold with conical curvature singularities at the world lines of the 
particles, such that these are timelike (lightlike) geodesies if the particle is massive 
(massless). It is possible to specify such a manifold as a geometric structure: M 
is locally modeled on M = M. 2,1 with transition functions in the isometry group 
ISO(2, 1) of M (M 2 ' 1 denotes Minkowski space, ISO(2,l) is the Poincare group). 
The holonomy map: 

p:7ri(M)->IS0(2,l) (1) 
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Figure 1: Examples ofag = 2,iV = 2 and g = 0, N = 4 surfaces E with the fat graph 
T above, the corresponding polygons P = E\r in the middle and the dual graphs 
in the bottom (Note that the dual graph of the genus two surface on the bottom 
left is reflected with respect to the plane of the paper). The few illustrated dashed 
lines are the intersections of the hyperplanes given by (J2J) with the t — spacelike 
planes. The edges of the polygon P are determined by the segments between the 
intersections of these lines. Those edges of the polygon bearing the same number 
have the same length and are to be glued with opposite orientation. The points A 
and B are examples of 3- vertices, Oj are 1-vertices. At vertices the sum of the angles 
are different from 2n. In the bottom figures the triangle bounded by the dual edges 
corresponding to those edges of V which are incident to the 3-vertex A are shaded. 
The dual graph cuts the surface into triangles and discs, corresponding to 3-, and 
1-vertices, respectively. 
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is associated to the construction: it is a homomorphism. We are considering globally 
hyperbolic manifolds with topologies M — E x I, for which 7i"i(M) = 7Ti(E), so the 
homotopy classes in M are identified with the homotopy classes in E. The holonomy 
map assigns Poincare group elements to each edge of 7. Since the number of its edges 
is more than the number of independent loops in a canonical homology basis of the 
fundamental group (E = 6g — 3 + 2N > 2g + N — 1), the set of generators is 
overcomplete, there are relations among them. 

Having the last paragraph said, we are ready to describe the explicit construction. 
Take a point Q in Minkowski space given by the origin of the Cartesian coordinate 
system X = (t,x,y). (It is a preimage of a point in M under the universal cover.) 
One of its images Vi(Q) = p(0)(Q) = Qi is given by the origin of another Lorentz 
frame X' = (ti,Xi,yi) such that X' = ViX. (Let us remind the reader that the 
dual edge Q is a closed curve, that is an element of the fundamental group.) Define 
hyperplanes in X by 

t = t i = {P i± Xf={Pt l xf. (2) 

These have spacelike normal vector. The intersection of these hyperplanes with the 
t — plane in X are straight lines. Choose a circular ordering of all sides of the edges 
of the fat graph V. The straight lines bound a polygon the edges of which are the 
intersections of the previously defined straight lines with their neigbours according 
to the circular ordering. For clarity let us summarize the procedure again. To each 
side of each edge of Y there is a corresponding oriented dual edge, to each oriented 
dual edge there is a corresponding Poincare holonomy, each Poincare holonomy maps 
point Q to point Qi and equations (J2J) and t = define a straight line. These straight 
lines bound a polygon 2 P which contains Q. See also edge 3 in the middle left part of 
Fig^for illustration. One can do the same procedure starting from Qi and obtain an 
isometric polygon meeting the previous one at edge i and having a timelike normal 
vector different from the purely timelike n = (1,0,0) of P containing Q. In fact 
each element of the fundamental group gives another copy of the polygon. This way 
we obtain a universal branched covering of the surface E with branch points at the 
vertices. For the case N = it is useful to visualize the situation as a nonplanar 
tiling of M 2 with the polygon P such that two copies meet at edges and three copies 
meet at vertices and the orientation is opposite for any two adjacent copies. The 
length Li± is the length of the edge given by the purely spacelike three-vector Ei± 
as the difference of the two vertices bounding the edge z± in coordinate system X. 
Li + = Li_ = Li by construction. The angles of the polygon can be determined from 
the consistency condition we have for each vertex of T. Assume that edge j, k, I are 
incident to a 3-vertex. Then 

PjPkPi = id, (3) 

2 Proving that starting from any geometric structure describing a globally hyperbolic Lorentzian 
spacetime with the given conical singularities one can construct a polygon without bad self inter- 
sections would amount to proving the statement that any universe with N ^ admits one 
polygon tessellation. We do not have such a proof. 
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since for the dual edges Q o ( fc o (; = e, where o denotes multiplication in the 
fundamental group and e is the unit element. Of course j above is either j + or j_ 
(and k, I likewise), one has to be careful about the orientation. After decomposition 
of the Poincare transformation as 

X, = V t X = AiX + at (4) 

with Aj G SO+ (2, 1) (+ denotes the identity component, the restriction means 
allowing only orientation preserving Lorentz transformations), and ^ 6 K 3 , we get 

Aj-AfcAj = id (5) 

and another relation, which involves also the a«'s. Using the definition of the boost 
and the rotation: 

/ cosh£ sinh£ \ / 1 \ 

5(0= 1 R(^)=\ cos^ -sin^ . (6) 

y sinh£ cosh£ / y sin^ cos ip / 

we can parametrize the Lorentz group elements as 

A = R^BprkWd (7) 

and the three angles at the vertex are given by 

a j = <Pi + 4>i a k = <P'i + 4>j ai = 4>'j + 4>k- (8) 

since the formula (j3J) rewritten with the help of (J2J) means that the parallel transport 
around a vertex can be written as a sequence of rotations and boosts. Note that 
the angle parameters <fi of the Lorentz group elements can be chosen without loss of 
generality in such a way that the angles a G [0, 2tt) at the vertices. The geometric 
picture is the following. The polygon containing Q meets the polygon containing Qi 
and that containing Qkoi at the vertex. If the normal vector of the first is denoted 
by n then A^n is that of the second and A^A^n = X^ l n is that of the third. The 
hyperbolic "angles" between the polygons at the edges are the boost parameters 
2r}j , 2r]k , 2rji. The relation n = A^A^A^n says that the parallel transport around the 
vertex does not change the vector n: space time is flat at the vertex. The group 
5*0+ (2,1) is the (orientation preserving) isometry group of the hyperbolic space 
H2 = {(£, x, y) : — t 2 + x 2 + y 2 = —1} with the induced metric from Minkowski space 
with constant negative curvature. The normal vectors are points in that space and 
the relation © give the angles tc — aj , tc — .tt — oi\ of a triangle 3 in terms of its 
lengths 2r]j ,2t]k ,27/;. We still need to specify the angles at 1- vertices. The normal 

3 The precise identifications of angles of a hyperbolic triangle in terms of the angles for the case 
when, say, aj > tt is explained in |25| . There can be at most one concave angle at a vertex [3J, 
which makes the functions a{rjj^rjk,rji) unambiguous for all three incident angles. 
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n is mapped to AjW ^ n under the holonomy, thus a 1-vertex carries nonzero (three 
dimensional) curvature. The angle there is 

& + #e[0,27r), (9) 

since the parallel transport around a 1-vertex is a boost and a rotation with the 
above angle. 

2.2 Constraints, degrees of freedom, symmetries 

In terms of the variables Aj and E^ we have constraints. 

• Half of the edge vectors are not independent since 

E i+ = -AiEi-, (10) 

where the minus sign is due to the opposite orientation. Using the decompo- 
sition of (J7J), one can write explicitly: 

= ±L i R{-4>' i )e^ E i+ = TLiR(<Pi)e {l) (11) 

with e W = (0,1,0). 

• There are Ag—2+N consistency conditions of the type of eqn.(jSJ) corresponding 
to 3- vertices. 

• There is a global constraint corresponding to the closure of the polygon P: 

+ £7«_) = 0. (12) 

i 

• There are N consistency constraints corresponding to 1-vertices. 

For the description of the latter we need the definition of the axis of a Lorentz 
transformation A. The axis is the eigenvector of A corresponding to eigenvalue 1 
and it reads: 

p = c(cosh 7] sin - - , — sinh^cos— — — , — sinhr/sin— — — ), (13) 
2 2 2 

and for later convenience we fix the constant c to be 

c — sgnsm — - — (14) 

whenever the rhs. is nonvanishing. If it is a timelike (lightlike) vector, the transfor- 
mation is called elliptic (parabolic). If it is elliptic, then it is conjugate to a pure 
rotation 

Ai = gR(2ir -m)g- 1 (15) 
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with g G 5*0+ (2, 1). The parameter m is interpreted as the mass 4 of the particle 
in units 8nG = 1, where G is Newton's constant. Taking the trace of (|15jh after 
elementary manipulations, we find: 

171 i a def / -t n\ 

cos — = ± cosh 7] cos (16) 

where a^e/ = 2n — (<p + <p') G (0, n) U (7r, 2tt) is the deficit angle at a 1-vertex with a 
corresponding elliptic Aj. The equation (jlfij) can be rewritten in the following way: 

p 2 + sin 2 ^ = 0, (17) 

where the square means the Minkowski scalar product with signature ( — h +). The 
above formula is the analog of the mass shell condition for the particle. To see that 
p is indeed the momentum of the particle, note that the spacelike direction of the 
velocity vector of the 1-vertex is given by 



E^-E i+ ~ (^ C os^—^^m^—^\ ~(0,p\p 2 ). (18) 



Taking into account eqn. ljlfij) we can conclude that the axis p is indeed the momen- 
tum of the particle 5 , the sign above gets fixed by choosing the orientation of the 
dual edge i. 

For a parabolic element, instead of formula (|T7jl . we have: 

p 2 = (19) 

this corresponds to a massless particle. However, the mass of this particle is not 
related to the deficit angle. On the other hand setting the parameter m = above 
does not corresponds to a massless particle, since it characterizes the trivial holon- 
omy A massless particle has a genuine parabolic holonomy, its deficit angle is given 
ClS db function of the boost parameter: 

cosh r] cos = ±1. (20) 

which is equivalent to the eq. ()19)) . The remaining N constraints are precisely the 
mass shell conditions and we assume that the Lorentz transformation A« correspond- 
ing to an edge incident to a one- vertex is never hyperbolic (which would be a tachyon 
with spacelike momentum). 

4 The definition of the mass comes from writing T 00 = m5(x — x(t)) for the static particle in 
the Einstein equations |22j . 

5 The unusual sine in the mass shell equation is a matter of convention It reflects the fact 
that 2+1 gravity is a toy model and the mass parameter coming from the stress energy tensor will 
appear with a unusual sine function in the mass shell condition. This analysis also shows, that 
arbitrary negative masses may be permitted since all the particles follow timelike geodesies. 
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Switching to the scalar variables L^r\x means that apart from the closure of the 
polygon P all other above enumerated constraints are solved identically. We use 
them to determine the redundant angle variables which determine the directions of 
the edge vectors as well. The closure of P remains as a constraint, it can now be 
written as 

^2LiZt(r,j) = 0, (21) 

i 

with the two vector Z{ = (Ei + + EiJ)/Li depending only on the boost parameters 
and 

= (2£-2)tt. (22) 

i 

restricting the boost parameters. In a different form of the latter after plugging in 
E = 6g — 3 + 2N one recognizes the Gauss-Bonnet theorem: 

H( Vj ) = 5> de/ = E ( E 27r - = 27r ( 2 " 2 9)- ( 23 ) 

v v \ieV J 

where the summation V runs in the set of vertices and i 6 V denotes the incidence of 
edge i to the vertex V. We have remaining symmetries corresponding to coordinate 
transformations of the chart: 



X i-> VX = AX + a (24) 

Vi i-> WiV 1 (25) 

The translation part of the Poincare transformation V has trivial effect on the 
variables {r]i, Li} and so does a rotation X i— > R(ip)X. So there are two independent 
nontrivial symmetries acting on the variables. One length parameter, say Li, can be 
chosen freely, it corresponds to a constant reparametrization of the time coordinate 
[2]. The number of physical degrees of freedom is the number of variables minus 
constraints minus symmetries: 

(2E- 1) -3-2 = 12g- 12 + 4iV, (26) 

the dimension of the reduced phase space, twice the dimension of the moduli space 
of iV times punctured surfaces 6 of genus g. 

The dynamics associated to the chosen time slicing can be determined from 
the equations (J2J). The boost parameters are associated to the hyperplanes defined 
by these equations, so they are time independent. The edge vectors E±± move 
perpendicular to themselves with a constant velocity tanh?7j. The variables are thus 
of action-angle type. We choose the convention that an edge i + (and its partner 
i_) moves into the polygon if rji < 0. Within a finite amount of time if an edge 
shrinks to zero length or a concave angle hits an opposite side, nine different kinds 

6 apart from the exceptions enumerated in the beginning in section 



10 



of transitions can occur 2\, the polygon may split into more polygons or a polygon 
may disappear, that is, the topology of the graph T may change. Accordingly, the 
number of variables may change as well. The values of the new parameters are 
determined unambiguously by the geometry, the new consistency constraints and 
the mass shell constraints. We shall not deal with this issue further in this paper 7 . 

Let us choose the Hamiltonian to be the total deficit angle given by (J23|) as 
a function of the boost parameters. This Hamiltonian generates the found time 
evolution [I] if the symplectic structure is 

{2r H ,L i ) = 5 i j. (27) 

3 First order formalism in 2+1 gravity 

In the first order formalism of gravity the physical fields are the triad e° and the spin 
connection a;" Locally both fields are one-forms with values in the Lie algebra of 
so(2, 1) = M 2,1 . Greek indices are tangent space indices, latin indices j, k... denote 
spatial components, latin indices from the beginning of the alphabet are indices 
from the Lie algebra, raised and lowered by the Minkowski metric with signature 
(—,+,+). To avoid confusion the index i will be used to index exclusively the 
finite number of variables corresponding to the edges of the graph V of the previous 
section. The (Hilbert-Palatini) action reads (with AttG = 1): 

I=\J ^ P Va b e«F b up , (28) 
where F stands for the curvature of the spin connection: 

F% = d »< - d »< + t abc u^„ c . (29) 

The theory does not exclude degenerate triads, hence it is more general than the 
second order formalism in terms of the metric 

g^u = VabC^l (30) 

However in case of nondegenerate triads the same equations of motion lead to the 
same solutions as those of the usual Einstein-Hilbert action and the symmetry struc- 
ture is also equivalent on shell. The physical interpretation for the triad is a local ob- 
server: in her Lorentz frame at the point x the spacetime vector v^(x) is e a ^{x)v^{x). 

In |25| it is proven for g > 1,N = universes, that there is an invariant smooth hyperbolic 
surface "triangulated" by 7 associated to each universe, which does not change during time trans- 
lations and symmetries. In that case there is only one kind of transition which can occur, and 
that corresponds to changing the triangulation by deleting one diagonal of a quadrilateral and 
drawing the other one. The new boost parameters and lengths are determined by trigonometry of 
the quadrilateral. Similar effect is expected to happen in the more general case N ^ at hand 

E2 
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With the help of the spin connection one can define the holonomy, that is the parallel 
transport operator along a curve a : [0, 1] — > M 3 : 

U a (a,b) =Vexp fj^ ds^^-cu^ , (31) 

where a(0) = a, a(l) = b, {T a ) b c = e b ac and V stands for path ordering. If v d is a 
Lorentz vector over b G M 3 then U a (a,b) c d v d is its parallel transported image over 
a G M 3 . In order to work in the canonical formalism, one decomposes the action a 
la Arnowitt Deser and Misner and arrives at the following form: 

S = JdtJ dFx^e^ie^ul + w$Dje b k + e a F b k ) (32) 

We assumed above the spacetime has the topology M = X x I where £ is a compact 
surface and / is an interval. The dynamical fields are the space components of the 
triad ef and the connection cu b , and one reads off from (jH2} that they are canonically 
conjugate: 

{e^x),cu b (y)} = e 0jkV ab 5(x,y), (33) 

where x and y denote coordinates on the surface S. The space of classical solutions 
is spanned by the the solutions of the constraints: 

F jk = 0, Dye h] = (34) 

which are equations of motions coming from the variation with respect to the La- 
grange multipliers eg and Uq, respectively. The torsion equations (Gauss constraint) 
generate the local 5*0(2, 1) transformations: 

uj(= oo a T a ) ^ g-Hg + g^ug, e a ^ (g' 1 )^ (35) 

and the curvature equations generate local translations: 

u i-> uj, e° ^ Dj\ a . (36) 

For a pedagogical account on the topic, see e.g. 

4 Reduction to covariant variables and adding par- 
ticles 

Our understanding of a point particle in 2+1 dimensions is based on the second order 
formalism. In [22] it has been shown, that with a stress energy of a massive point 
particle the solution to the Einstein equations is a cone. In Cartesian coordinates 
a particle is a point in Minkowski space with a wedge cut out and its edges glued 
together with time translation proportional to the spin of the particle and the angle 
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of the wedge is proportional to the mass. To incorporate point particles in the 
theory the first step is to add the necessary terms to the action (J32J), such that we 
obtain the correct equations of motion for them. The analysis has been done in 
both the mathematics and physics literature [23 IH1 EH EE EEJ • Below we briefly 
summarize the procedure and the reduction of the action to a finite number of 
covariant variables. The references f3J E] are followed with some modifications in 
order to incorporate also higher genus surfaces as spacelike slices, that is, spaces with 
nontrivial tangent bundle. The last step, which is the main result of this paper will 
be done in section arriving at the symplectic structure of the polygon model and 
its Hamiltonian as the constraint (s) remaining from the reduction. If there are point 
particles present, the first order and second order gravity are not even classically 
equivalent jjj, due to the different structure of large gauge transformations versus 
large diffeomorphisms. We do not deal with this fact here. We also restrict ourselves 
to the case of spinless particles. 

Consider the equal time surface E as a simply connected region P = S/T, with 
the graph T of sectional and pairwise identifications of the boundaries (which are 
the worldlines of the edges connecting vertices of T). On shell, the spin connection 
and the triad are given by pure gauge [Ti] : 

uj« b = (g-%gy b e» = (ff- 1 )W (37) 

with potentials g : R 2 x / -»• SO (2, 1) and / : R 2 x / -> R 2 * 1 . They are defined on 
the universal cover, they are not single valued on the surface S. However, the fields 
should be single valued on £ which restricts the potentials. That is, most generally, 
on each pairs of edges i + and z_ in dP which correspond to the same edge i of V 
the following condition should hold: 

g l+ = A l9l ., fi + = ^ibU'i- ~ a \) (38) 

with Ai G 50(2, 1) and a { e R 2 ' 1 - We need to impose the consistency condition: 

AiAjAjt - id = (39) 

for all triples k) corresponding to 3- vertices of T. This way the constraints in 
(jH2j) are solved everywhere, but at the 1-vertices. We want the 1-vertices to be point 
particles. To this end, we have to impose extra constraints: 

trAj - 2 cos rrij - 1 = (40) 

where i runs in the set of labels of 1-vertices. They make sure that the massive parti- 
cles move along timelike geodesies, massless particles move along lightlike geodesies 
and there are conical singularities at their locations ^3JE3- A second set of variables 
is given by the following formula: 

E? ± = [ C±- (41) 

Ji± 
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The potential / defines an embedding of £ to Minkowski space, and the "edge 
vector" Ef ± is the relative position vector in the background Minkowski space of the 
two ends of edge i± of T. The image is a surface bounded by a (nonplanar) polygon, 
with pairwise identifications of the edges. The edge vectors are not independent, 
but satisfy 

E? + = -A? b El, (42) 

where the minus sign is due to the opposite orientation of the gluing. Finally there 
is a global constraint they also obey: 

i 

The computation of the symplectic potential is now identical to page 66 of [To] . 
Instead of repeating the calculation, we enumerate the notational differences. First, 
we have only one polygon, so g/\ is just the potential g. Then g\ <-> Aj, Matschull 
labels the edges of V with A and z\ <-> Ei- with A > 0. Finally we use the vectorial 
representation of SO (2, 1), while he uses the spinorial. The result of the reduction 
reads: 

9 = I dV Vab dw?e> = -J2 <L{K l dK)\El. (44) 

The above equation is the result of a page of calculation: Matschull plugs in the 
expressions (pTTjl. integrates over P, but the integrand is a total derivative. Hence, 
the two-dimensional integral becomes a sum over integrals for edges, which is again 
a total derivative and one finally obtains the rhs. of ([44)1 . We will need the explicit 
expression for the Poisson bracket of the components of the global constraint ([43)1 . 
Since only the following brackets are nonzero among the edge vector components 

m 

{E? ±1 El} = e?E^ (45) 

the global constraints (|4*H|) generate Lorentz transformations of the edge vectors and 
obey the same algebra: 

{C\ C b } = ef C c . (46) 

The brackets of the rest of the constraints are analyzed in the literature ^SIE], it 
will not be repeated here. Instead we go on with the symplectic reduction to recover 
the symplectic structure of the polygon model. 



5 't Hooft polygon model in the first order for- 
malism 

Waelbroeck and Zapata have shown j2Ij that a triangular lattice model based on the 
first order formalism also reduces to the polygon model if one introduces the scalar 
variables of the polygon model from the covariant parallel transport matrices and 
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edge vectors, which solves some constraints identically. We will now show how the 
't Hooft polygon model arises in the smooth first order formalism. First we require 
by means of a gauge transformation that the initial time surface is mapped to a 
planar polygon by the potential / to Minkowski space. Then the scalar variables 
will be introduced by means of solving the constraints introduced in the previous 
section. One can proceed in two ways to recover the symplectic structure (J27)l of 
the scalar variables of the polygon model. The first is further reduction of the sym- 
plectic potential given by formula (jH)) . The second is the derivation of the Poisson 
brackets directly by appropriately expressing the scalar variables as functionals of 
the spin connection and the triad. Finally in the last section we identify the induced 
constraints, which provide dynamics for the model. 

5.1 Phase space 

The edge vectors Ef± and the holonomies A, are the finite number of covariant 
variables which contain the degrees of freedom of the system. To get to the polygon 
model we perform a local translation of (jBEJ) to obtain a planar polygon. Each edge 
vector Ef± has to lie in a spacelike hyperplane in Minkowski space. Denote the unit 
timelike normal to this hyperplane by n a . Let us define the boost parameters 8 and 
the lengths: 

cosh277 i = n A> 6 . (47) 

U = yjE? + El Vab = y/Ef_El Vab (48) 

Now we can compute the symplectic structure departing from the rhs. of (|44p. Let 
us choose a coordinate system in which the normal n a is purely timelike. The edge 
vectors are purely spacelike and we can use the explicit form given by (fTTjl if A, is 
written as in eqn.(JZJ). Plugging this into the term e b ac (hj l dAi) a b Ef_ the result of a 
direct calculation shows that the terms dip and dip' vanish and the coefficient of the 
d(2r)i) term turns out to be —Li. Thus the symplectic potential is as expected: 

e = -2j2L i dr li (49) 

i 

To understand better the origin of the result above, we present an alternative deriva- 
tion of the Poisson brackets directly from the original fields of the theory. The 
Lorentz group element A» is the holonomy along the dual edge Q. More precisely, if 
the points a_ and a + are the starting and endpoints of a curve freely homotopic to 
the dual edge £ with the same orientation, then the following relations are true: 

[/ ( a 6 (a_, a+) = g-\a_)g(a+) = g- 1 {aJ)Kg{aJ) = g- 1 (a+)Ag(a+). (50) 

8 Let us remark that for an edge i, it is possible to choose the gauge in such a way that Ai is a 
pure boost hence trAj = 2 cosh 2r)i + 1 as proposed in [2Hj , but this cannot be done globally for all 
edges. It would mean that the (cosh of) the boost parameters are equal to traces of holonomies 
along closed curved, thus gauge invariants. We will see that this is not the case 
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Note that the index i is omitted from above to keep the notation simple. The 
holonomy can be also expressed in the usual way as the functional of the spin 
connection given by eqn.([3ip. From now on we choose a fixed edge i of V and a dual 
edge of 7 denoted by £. The expression for the Poisson bracket of the associated 
length and boost parameter reads: 

{L, cosh 277} = {L,n a (a + )U* b (a_,a + )n b (a + )} = 

e i,n cd n (a )n»(a ) f £x 5Li 6U ^ a ~^ (51) 
e 0jk V n a {a + )n (a + ) J^d x ^ ^ 

where n a (a + ) = g~ 1 (a + ) a b n b . In general, the functional derivative of the holonomy 
can be written as follows: 

5U^ f , dxHs 



r dxHs) 

where (1 is the segment of ( until the point x, (2 is the segment of ( from the point 
x and x(s) is the parametrization of the curve (. The variation of the the length 
with respect to triad can be written using equations (|3T|). (jUJ, as 

5L Eb + f b dx^ 



5ej(x) 



drg(r) b c ^-5(x(r),x), (53) 



where x{t) is a parametrization of the edge. Now we substitute (jB^j) and (jo^j) to 
(j51|) and integrate out one Dirac delta. If the integral for ( is done first from the 
point S- = x{t) to its image s + along the curve homotopic to (, then we have (1 = 
and (2 = C> so we can write: 
{L, cosh 2r/} 

dx j (r) dx k (s) 



dr / ds 



C 



Eb+9 [ S+)b ° " n c (s + ) (T a U(s-iS + )) c b n\s + ) 



(54) 

The expression outside the square bracket is independent of s, since we can rewrite 
it using the notation {v x w) a = e a b c v b w c as 

g-^s^E+ng-^s^n x g~\s + )An) a = E a + (n x An) a , (55) 

because both the vector product and the Minkowski scalar product are invariant 
under Lorentz transformations. The integral of the square bracket in equation (j54j) 
is the homotopy invariant oriented intersection number of the edge i and the closed 
curve (. It is one for the pairs (i, Q) and zero otherwise. The last ingredient is to 
verify the following formula: 

n x AiU = smh(2i]i) Ei + . (56) 

E i+ is clearly orthogonal to both n and A^n due to (pf2*|) and the scalar factor in the 
previous equation is also not difficult to check using the definition (|4Tj). We have 
thus arrived at the desired result: 

{Li, cosh 2r]j} = 5ij sinh 2^ — > {Li, 2r)j} = Sij (57) 
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5.2 Dynamics 



At the intermediate stage of the symplectic reduction, the finite number of covariant 
variables Aj and Ei have to satisfy a number of constraints, given by (J39|) . pH|l . (jH)) 
and (|4Hjl . Changing to the scalar variables Li , r]i identically solves the first three, if 
the triangle inequalities 

|»7i| + \Vj\ > \Vk\ (58) 

are satisfied, for every vertex V = {i, j, k}, and for every permutation of {i, j, k}. The 
explanation is the following: Eqn. (j39|) and (j4*Uj) determine the angles of the polygon 
in terms of T]i, eqn. ()42j1 determines the direction of the edges. In reference j2I] it is 
argued, that the Gauss constraint is solved by introducing the scalar variables and 
there remain induced curvature constraints, but the situation is on the contrary: 
the Gauss law and only that remains after the reduction. Namely, the closure of 
the polygon P: (J4*3j) written as (|2ip with only two independent components and the 
sum of angles in the Euclidean polygon equivalent to the Gauss-Bonnet theorem for 
the surface it represents, given by eqn. (|23jl . The latter is the Hamiltonian, which 
generates a time evolution that preserves the closure of the polygon by construction 
H|, thus 

{H,C a } = J2 L i z i = ^ (59) 

i 

where the dot now indicates the time evolution induced by H. The fact that 
{C a , C b } = e ab C c does not get modified is difficult to check explicitly. However, 
the correct number of degrees of freedom comes out as in eqn. (|2fij) only if all three 
independent constraints H, C\ and Ci are first class and generate time translation 
and two independent Lorentz transformations respectively (since a rotation is fac- 
tored out by using lengths and angles instead of edge vectors). In other words the 
constraint algebra closes and the most general time evolution is generated by the 
following Hamiltonian: 

H' = aH + bC l + cC 2 . (60) 

6 Discussion 

In this paper, the 't Hooft model of globally hyperbolic 2+1 dimensional Lorentzian 
universes with zero cosmological constant and compact spacelike part has been de- 
rived in terms of the fields of the first order formalism of gravity. We could reduce 
the originally infinite dimensional phase space spanned by the spin connection and 
the triad to the space of finite number of scalar variables in a simple way for the 
generic case of arbitrary genus and number of point particles. The induced symplec- 
tic structure has been derived both directly from the smooth fields and from a finite 
set of covariant variables of Matschull, Waelbroeck and Zapata without the need of 
introducing any lattice. Apart from showing the connection of different approaches 
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to 2+1 gravity, interpreting the polygon model directly as a gauge theory has an- 
other advantage. While the original treatment in terms of geometric structures has 
a "rigid" time evolution generated by H and induced by the slicing of space time, in 
gauge theory, there is an arbitrariness of taking any linear combination of the con- 
straints to be the Hamiltonian. One may use this freedom to prove that one-polygon 
tessellation is generic and transitions corresponding to splitting (and disappearance) 
of polygons may be completely avoided. Another aspect which makes it worths to 
pursue progress in this direction is the cotangent bundle structure of the polygon 
model. It is a convenient starting point for a quantum theory. However this has only 
been shown and an explicit section given for the case of no point particles. 
For the case of point particles there is no complete classification of the phase space 
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